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ABSTRACT 

As the next step toward an improved large scale Galactic magnetic field model, we 
present a simple comparison of polarised synchrotron and thermal dust emission on 
the Galactic plane. We find that the field configuration in our previous model that 
reproduces the polarised synchrotron is not compatible with the WMAP 94 GHz 
polarised emission data. In particular, the high degree of dust polarisation in the 
outer Galaxy (90° < £ < 270°) implies that the fields in the dust-emitting regions 
are more ordered than the average of synchrotron-emitting regions. This new dust 
information allows us to constrain the spatial mixing of the coherent and random 
magnetic field components in the outer Galaxy. The inner Galaxy differs in polarisation 
degree and apparently requires a more complicated scenario than our current model. In 
the scenario that each interstellar component (including fields and now dust) follows 
a spiral arm modulation, as observed in external galaxies, the changing degree of 
ordering of the fields in dust-emitting regions may imply that the dust arms and 
the field component arms are shifted as a varying function of Galacto-centric radius. 
We discuss the implications for how the spiral arm compression affects the various 
components of the magnetised interstellar medium but conclude that improved data 
such as that expected from the Planck satellite will be required for a thorough analysis. 

Key words: ISM: magnetic fields ~ Galaxy: structure - polarisation - radiation 
mechanisms: general - radio continuum: ISM 



1 INTRODUCTION 

Determining the large-scale structure of the Galactic mag- 
netic fields is a longstanding and difficult task. The fields 
are important to a variety of astrophysical phenomena from 
star formation to cosmic ray propagation, but there are few 
ways of observing them directly. Radio astronomy is the pri- 
mary way that we study magnetic fields, but any analysis 
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in only one domain is subject to a variety of degeneracies 
that cannot be disentangled without several complementary 
measurements and a multi- wavelength approach. 

IJaffe et all l|2010l . 120111 ) model the large-scale Galac- 
tic magnetic fields in the plane using a combination of 
three observables: the total and polarised intensity of dif- 
fuse synchrotron emission, and Faraday rotation measures 
(RMs) of extra-galactic sources. Other analyses of the 
large-scale ma gnetic field s have been performed us i ng the 
full sky (e.g ISuii et al.1 |2008|. iFauvet et al.l 1201 ll . |2012| . 
Ijansson et al.J 200E: Ijansson fc Farrarl 20121 ) or the high- 
latitude sky |strong et al.ll201ll ). We focus on the plane in 
order to probe through the full Galactic disk in detail, and 
we now continue these studies with the addition of the po- 
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larised thermal dust emission. In addition to being an in- 
dependent tracer of the field direction, the dust emission 
follows a different spatial distribution and therefore is not 
subject to the same uncertainties as previous studies using 
only synchrotron emission. 

An important aspect of our previous work is to sep- 
arate the fields into components defined by their observa- 
tional signatures. As described in I Jaffe et al.l (|201Cll ). the co- 
herent fields contribute to both polarised emission as well as 
to RM, while isotropic random fields contribute only to the 
total intensity of the diffuse emission. A third component 
is needed to include the effect of anisotropy in the turbu- 
lence, a component that maintains a common ordered axis 
but changes direction stochastically. This "ordered random" 
field component contributes to polarised emission but not to 
RM.0 

Since the cosmic ray electron (CRE) distribution is 
thought to be fairly smooth (see §3] for a discussion), the 
contributions to the synchrotron emission from the different 
magnetic field components do not depend strongly on their 
relative spatial locations. Our aim is therefore to add the 
information from the polarised dust emission from WMAP 
to study the spatial distribution of each component. The ge- 
ometry of how the polarised emission relates to the different 
magnetic field components is identical for diffuse thermal 
dust emission. But the distribution of dust is not as smooth 
as the C RE distribution and is know n to peak in spiral arms 
(see, e.g.. lDrimmel fc S pcrgcl 200IV This allows us to probe 
whether the field components peak i n the same ridges or 
whether, as seen in external galaxies l|Beckll2009l ). different 
components may lie along different arms that may have dif- 
ferent pitch angles or be shifted relative to each other. Note 
that WMAP W-band at 94 GIfz primarily probes the rel- 
atively cold dust, which has a different spatial distribution 
from the hotter dust dominating the infra-red. 

We hope also to shed light on the outstanding question 
of arm versus inter-arm contrasts in the components of the 
magnetised interstell ar medium (ISM) . In external galaxies, 
iFletcher et al.l (l201ll ) and references therein do not find the 
contrast in synchrotron emission between arms and inter- 
arm regions as predicted by the scenario of compression in 
the spiral arm shocks. It remains unclear if this is due to 
the limited resolution of the current observations in external 
galaxies. Our method comparing the complementary observ- 
ables of synchrotron, dust, and RMs along with a variety of 
external constraints allows us to ask similar questions from 
within the plane of the Milky Way. 

We test a model of the magnetised ISM where each com- 
ponent contains a spiral arm structure but the arm peaks 
are not necessarily spatially coincident. In particular, we ad- 
dress the question of how to reproduce the level of dust po- 
larisation in the outer Galaxy (defined as 90° < £ < 270°, 
i.e. quadrants two and three). This region is dominated by 



^ This component is rarely treated by itself, and there is some 
confusion in the li t eratu re over how to unambiguously refer to it. 
Ijansson fc Farrail | |2012| '| refer to it as the "stria ted" compo nent. 
For consistency with previous literature (e.g.. iBeckl 120091 ') and 
to remove any ambiguity, we prefer the term "ordered random", 
noting that simply referring to the ordered field, as we have done 
in the past, might also be mistakenly construed to include the 
coherent component. 



the Perseus arm that passes outward of the Sun's position. 
The outer Galaxy is therefore ideal to model unambigu- 
ously the arms of the different components, since we are 
looking through one inter-arm region and one arm ridge, 
beyond which is expected to be very little ISM. The in- 
ner Galaxy (defined as —90° < £ < 90°, i.e. quadrants one 
and four) is more complicated, and there are already indi- 
cations th^t_tlMj^xge;;Scale field structure is not a simple 
spiral (|van Eck et al.ll201ll ). 

Our simple comparison of polarised synchrotron and 
dust emission in the outer Galaxy allows us to study the 
relative positions of the arms as well as the arm versus 
inter-arm contrasts. T his is analogous t o stu dies such as 
IPatrikeev et al.l (|2006l ) or iFletcher et~ai1 (|201ll ) perform on 
M51 but with a very different perspective from within the 
disk. This is important to understand the relationship be- 
tween the magnetic fields and other components of the ISM 
and the effect of spiral arm compression on all components. 
We can test the scenario where gas, dust, and fields are 
compressed in the spiral arm shock, which predicts strongly 
ordered fields in the dust-emitting regions. We further dis- 
cuss the implications of these results for the inner Galaxy 
and outline future work. 

We describe the data we use in § [5] and the modelling 
and simulation methods in § [3] We then discuss the inputs 
and existing constraints on the various components of the 
magnetised ISM in § [J] We show the results of our simula- 
tions in § [5] and discuss the implications in § (6] along with 
outlining the further work these results point toward. 



2 DATA 

Following our previous work, we use a combination of low- 
frequency radio as well as microwave data, in total and po- 
larised intensity, and examine the profiles in longitude of 
each frequency along the Galactic plane. 



2.1 Radio and low- frequency microwave datasets 

The datasets and referenc es are summar i sed in Table [1] We 
use the 408 MHz survey bv lHaslain et all (1 19821) . the W MAP 
23 GHz polarisation map from Jarosik et al.l (l201ll ). and 
three catalogues of RMs, the Southern a nd Canadian Galac- 
tic Plane Surveys (SGPS and CGPS) of lBrown etall (|2003| 
2007!) and the additional catalogue from Ivan Eck et al.l 
(|201ll ). 

Each of the emission maps in HEALPijfl i Gorski et al.l 

(2005!) format is smoothed to a common low resolution of 
FWHM 6° and downgraded to iVsido = 32. The RM cat- 
alogues are averaged into bins at the same resolution and 
assigned error bars based on the RMS variation of the in- 
dividual measures within the bin. The profiles in longitude 
along the plane are shown in Fig. [T] along with a comparison 
to an updated form of our previous model (| Jaffe et al.ll201lh 
as described in § 13.11 (Note that we do not attempt to fit 
the region within 10° of the Galactic centre, and the model 
is not plotted there.) 

In order to use the low-frequency radio observations in 
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total intensity on the plane, we apply a correctio n to remove 
most of the free-free emission. As described in Ijaffe et"al] 
(j201ll ). we find that the WMAP MEM free- free solution at 
33 GHz appears to overestimate the thermal contribution. 
We found, using simple cross-correlations, that a normaliza- 
tion factor of 0.8 appears to give a good correction to the 
radio frequencies, and we use it in the dust band total in- 
tensity at 94 GHz as well. In both cases, the template is 
extrapol ated to each corr e cted f requency with the spectrum 
given in iDickinson et al. I (|2003'). This component has been 
subtracted from the total intensity profiles in green in Fig. [T] 
and is also shown in orange for comparison. 

Note that the polarised intensity, PI = \/Q^ + plot- 
ted in Fig. [T] is subject to a noise bias due to the squaring 
of the noise on Q and U. In the case of the synchrotron 
polarisation on the plane, the polarisation bias is too small 
to see. But for the dust polarisation at 94 GHz, this bias is 
significant. To reduce it, we combine the data from different 
WMAP years into two independent maps each of Q and U, 
and then construct xPI = \/QiQ2 -I- UiUi- Since the noise 
contributions in Qi and Q2 are independent, this cross es- 
timator of the polarised intensity has little noise bias and 
roughly equivalent uncertainties. A better estimator of PI 
would depend on the detailed properties of the detectors, 
and for the quantitative analyses in our previous work we 
used Q and U separately, but for the qualitative analysis 
presented here, this xPI is quite informative. 



2.2 WMAP 94 GHz polarisation 

The new dataset we include in this work is the WMAP 
94 GHz frequency that, in polarisation, is dominated by 
thermal dust emission. As can be seen in Fig.[T] the free-free 
emission correction is also necessary on the plane, but unlike 
at 23 GHz, it does not dominate the total intensity. 

The calibration and signal-to-noise ratio of the 94 GHz 
polarisation is more of a concern, however, since the dust 
emission is more weakly polarised than synchrotron emis- 
sion. As noted above, the noise bias is significant at 94 GHz, 
which we overcome by plotting xPI. We can also charac- 
terise the noise and calibration uncertainties by comparison 
of the maps for individual years. We use the variance from 
year to year in Q and U to estimate the uncertainty on xP7, 
which then includes the effects of both the noise and the sys- 
tematic uncertainties that are of that time-scale. These un- 
certainties are plotted as the dotted green lin es around the 
green xP/profile in Fig.lTl ljarosik et al.l |20l3) single out the 
W4 difference assembly as particularly prone to systematic 
artefacts, and we find that including it does significantly in- 
crease our estimate of the uncertainties. We therefore follow 
Jarosik et al. and omit W4 from the analysis. 

The S/N apparent in Fig. [T] is somewhat low for per- 
forming quantitative fits and parameter exploration but is 
sufficient for the simpler qualitative analysis we present here. 
For a thorough exploration of the parameter space, we will 
need be tter data such as that expe cted from the Pianc/Jf] 
mission ijPlanck Collaboration! I2OIII ). 



3 METHODS 

T he methods and t o ols used in this work are described 
in iJaffe et all (|201(]| '). Ijaffe et al.l l|201ll) , and iFauvet et all 
In short, we have adopted models for each of the 
physical components of the ISM (thermal electrons, mag- 
netic fields, dust, and CREs) and then integrated through 
the Galactic plane to simulate observable polarised emission 
and Faraday RM for comparison to the data. 

As discussed in § [2]2l the WMAP 94 GHz polarisation 
data do not have a high enough S/N for a quantitative anal- 
ysis such as the MCMC parameter exploration. So for this 
work, we simply compare the profiles and draw qualitative 
conclusions from the morphology. 



3.1 Models of the Galaxy components 

We use a magnetic field model based on a spiral arm com- 
pression as described in jjaffe et al.1 (^010). We have mi- 
proved this model somewhat to simplify the parametrization 
and to make the field continuous and closer to divergence- 
free. See Appendix |X] for the modified parametrization. The 
principle differences are: to remove the extra complexity of 
the radial profile in the innermost few k pc of the Galaxy 
that orig inally followed iBroadbentI l| 19891 ) but is excluded 
in this analysis; to define the transition between arms con- 
tinuously (rather than to jump discontinuously in the inter- 
arm region as in the previous parametrization); to vary the 
width of each arm with Galacto-centric radius in inverse pro- 
portion to the coherent field strength. Remember that the 
coherent field is assumed parallel to the arm, so that the 
arm corresponds to a magnetic fiux tube. This last modifi- 
cation, which conserves magnetic flux, maintains a roughly 
divergence-free field. We also allow the spiral arm patterns 
for the different components (which are otherwise similar) to 
vary independently in orientation about the Galactic north- 
south pole. In other words, the components' arms can be 
shifted rel ative to each othe r by a constant angle. 

As in jjaffe et all l|201lh . we use the "NE2 01" t hermal 
electron distribution model of ICordes fc Laziol 1I2OO2I ). 

We also use the GALPROi|j cosmic ray propagati on code 
with a model based on the "z04LMPDS" from Strong et al.l 
(l20ld) Fl a nd w ith the fitted injection spectrum from 
iJaffe et al.l (|201ll ). Based on the given distribution of CRE 
sources, the code determines the spatial and spectral dis- 
tribution of CREs after propagation in the Galaxy's mag- 
netic fields, in this case our model rather than the GALPROP 
model. The resulting distribution then is taken to compute 
the synchrotron emission as described in § 13.21 

We use a dust distribution model based on the same spi- 
ral arm parametrization used for the magnetic fields, which 



|http: //planck. esa. Int I 
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5 The Z04LMPDS was modified in Jaffe et al. (2011) to reduce the 
grid resolution for speed, to disable the irrelevant nucleon propa- 
gations, and to change the magnetic field used for the synchrotron 
losses. Most of the p arameters for z04LMPDS are defined in ta- 
ble 1 of lStrong et al.l ((2010 j with t he exception of the C R source 
distribution. The formula is that of Lorimer et al. 1 1I2OO6I ') (Eq. 10) 
but adjusted in IStrong et al.l 1I2OII ) to agree better with Fermi- 
LAT gamma rays, namely B = 0.4751 and C = 2.1657 and con- 
stant for R=10-15 kpc. 
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Figure 1. Profiles of tiie observables {thick green lines) and tiie sligfitly updated version of our previously fitted model {thin blue lines) 
along the Galactic plane. See text in § [21 In order from the top: total intensity at 408 MHz minus free-free estimate; polarised intensity 
at 23 GHz; model synchrotron polarisation fraction; total intensity at 94 GHz minus free-free estimate; polarised intensity at 94 GHz 
and (as the pale green filled region) its uncertainty (see text in S I2.2|l : dust polarisation fraction; Faraday RM at 1.4 GHz. The shaded 
regions show the Galactic variance for the model, i.e. the expected variation due to the turbulence. The dashed line in the dust panels 
shows the modelled synchrotron contribution at 94 GHz, which is small but not negligible. The dot-dashed gold line shows the free-free 
estimate compared to the total intensity curves. 



is in fact similar to that of iDrimmel fc Spergell (|200ll ). Our 
parameters, given in Table 1X21 are adjusted to match the in- 
tegrated total intensity profile of the data at 94 GHz. Fig.[T] 
shows that the model profile matches well at large scales. 

We use the dust polarisation model described in 
iFauvet et all l|201lf) . i.e. with a geometrical factor of sin(Q) , 



where a is the angle between the local field orientation and 
the line-of-sight (LOS). This factor approximates the effects 
of grain alignment and projection onto the LOS, but we note 
that our analysis is not particularly sensitive to this factor. 
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Survey 



Frequency [GHz] Original resolution 



Coverage 



Reference 



Haslam 
WMAP 7-year 
WMAP 7-year 
CGPS RMs 
SOPS RMs 
VLA RMs 



0.408 
22.8 
94 
1.4 
1.4 
1.4 



~51 arcmin Full sky 

53 arcmin Full Sky 

13 arcmin Full Sky 

1 arcmin part of Northern plane; see Fig. [T] 

100 arcsec part of Southern plane; see Fig. [T] 

50 arcsec part of plane; see Fig. [T] 



Haslam et aL (1982) 
Jarosik ct ak (2011) 
Jarosik et al. (2011) 
Brown et al. (2003'! 
j3rown et al. (2007) 
va n Eck et al.^ (2011) 



Table 1. The datasets used in this analysis. The RMs are the extra-galactic sources only from the Canadian Galactic Plane survey 
(CGPS), the Southern Galactic Plane Survey (SGPS) and the van Eck et al. survey using the Very Large Array (VLA). 



3.2 Simulation of observables 

As described in IJaffe et al.l ()2010l l. we simulate the diffuse 
polarised emission a t large scales using the publicly available 
code HAMMURAslfl (jWaelkens et al.l [20091 ). which simulates 
observables such as synchrotron and dust emiss i on in full 
Stokes parameters and Faraday RMs. Ijaffe et al.1 (|201ll ) de- 
scribed how we then linked this code to the Galprop cos- 
mic ray propagation code in order to include high-energy 
constraints from both direct cosmic-ray measurements or 
from diffuse 7-ray emiss i on as measured b y the Fermi LAT 
l|Ackermann et al.ll20ldl : lAbdo et al]|2009l ). This is impor- 
tant to compare the synchrotron total intensity at low radio 
frequencies with polarised intensity at microwave frequen- 
cies. (Since low frequency polarisation is affected by Faraday 
depolarisation and higher frequencies by free-free contami- 
nation on the plane, it is complicated to compare total and 
polarised intensities on the plane at the same frequency.) 
Parameters have simply been adjusted manually to match 
the data as plotted in profiles along the plane as in Fig. [1] 
(A quantitative parameter exploration will be performed in 
a future work.) 

Fauvet et al. use HAMMURABI to compute the "geo- 
metric suppression" map, i.e. the map of the polarisation 
fraction, and multiply it by a total intensity template to 
scale the polarisation simulation. This allows a more ac- 
curate modelling, since the intensity template partly ac- 
counts for localised features not in the model. In Fig. [2l 
we follow the same approach, using the 94 GHz total inten- 
sity data correcte d for the free-free emission as described in 
iJaffe et al.l (|201ll ). This is in contrast to Fig. [T] where we do 
not use a template but simply plot the result of the integra- 
tion through our dust density model. In this case, both the 
data and the simulated total intensity profiles are plotted for 
comparison, and their very similar large-scale characteristics 
give confidence that the dust distribution model should be 
accurate enough for this analysis. 

As described in IWaelkens et~ai] (120091 ') . the simulation 
resolution in HAMMURABI is refined with radius as the ob- 
served solid angle encompasses a larger volume. The result 
is a resolution that varies but does not increase significantly 
with radius. In our large-scale, low-resolution simulations, 
it is roughly 20 pc or larger. With an outer scale of turbu- 
lence of 100 pc, this is clearly not sensitive to more than the 
largest scales of the turbulence (represented by a Gaussian 
random field simulation). Any structure on scales smaller 
than this will be modeled only as an average over the "sub- 
grid" effects. We have tested that increasing the resolution 

^ |http : //sourcef orge .net/projects/hammurabicode/l 



does not impact our results. The maps themselves are simu- 
lated to match observations smoothed to 6° as described in 
§ [2] above. 



4 INPUTS AND CONSTRAINTS 

We use our previously studied models and additional infor- 
mation from other studies that is relevant to the comparison 
of data and model. In the following, we explain the choices 
we have made and the relevant external constraints. 

4.1 Synchrotron 

In IJaffe et al.1 l|2010l , HoTl), we constrained the ratios of 
the magnetic field components (defined above) using syn- 
chrotron total and polarised emission as well as Faraday 
rotation measure. Fo r the curren t analysis, we retain largely 
the same model as in IJaffe et al.l (|201lh but allow for an az- 
imuthal shift between the spiral patterns of the different 
components. This has little effect on the integrated syn- 
chrotron profiles, since the model CRE distribution is rela- 
tively smooth and varies slowly with Galacto-centric radius. 

4.2 Arm versus inter-arm contrast 

The theory of shock compression at a spiral arm predicts a 
gas density contrast of approximately a factor of 4 in the 
shocked versus the pre-shock medium, assum ing a strong 
adiabatic shock (e.g.. iLandau fc Lifshit j|l987l ). The compo- 
nents of the magnetic field that are parallel to the shock are 
therefore compressed as well and thus amplified by a factor 
of 4. This is the physical motivation for the creation of the 
ordered random field from the isotropic random field and for 
the amplification of the coherent field, which is parallel to 
the same spiral shock front in our model. 

We do not model the shock itself but simple spiral 
arms with Gaussian profiles. This can be thought of as 
a large-scale approximation to the modulation of a shock 
that then re-expands to its pre-shock strength downstream. 
We can define a contrast as the ratio of the peak in the 
arm to th e minimum b etween arms. The original model of 
iBroadbc nt et al.l (|l990[ ) used a contrast of 3.5 for both co- 
herent and isotropic random field components, which they 
found reproduced well the synchrotron total intensity pro- 
file, particularly the amplitude of the emission observed in 
the direction of spiral arm tangents. 

We use the Broadbent et al. value of 3.5 as a starting 
value, but in our model, the strength of the co herent field 
component in each arm is allowed to vary as in IJaff'e et al.l 
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l|2010l ) so as to be consistent with the RM data. This con- 
trast then varies for each arm by the factor ai given in 
Table lAll But these values cannot be compared with the 
theoretical value of 4 for shocked versus pre-shock medium 
without some knowledge of the width of the shocked region 
and its re-expansion, since our Gaussian profiles are effec- 
tively a large-scale average. Note also that the coherent field 
strengths in the arms we choose to match the RM data are 
degenerate with the thermal electron distribution, which is 
not very precisely known. 

The ordered random component increases at the shock 
as a result of compression of both any pre-shock ordered 
random field and of the pre-shock isotropic random field. It 
then decreases back to its pre-shock strength, due to both 
re-expansion across the spiral arm and re-isotropisation by 
supernova-driven turbulence after star formation triggered 
by the shock. There are also other possible sources of the 
ordered random component such as superbubbles or shear- 
ing motions. In our model, this component has a Gaussian 
profile that goes to zero between the arms (and therefore 
has no meaningful contrast), so this is clearly only a simple 
approximation to include a variety of physical effects. 

The isotropic random component itself would remain 
unchanged through the shock, only providing the source of 
an ordered random field created there. Our modeled spiral 
arm modulation of the isotropic random component there- 
fore has a difi'erent origin such as the supernova-driven tur- 
bulence downstream of the shock. The arm contrast is there- 
fore not constrained a priori and is left at the Broadbent et 
al. value of 3.5. 

The diffuse component of the gas and dust would also 
be compressed by a factor of 4 in the shock, but the aver- 
age gas and dust density is more strongly compressed due 
to the formation of dense clouds behind the shock. Since 
the WMAP data at 94 GHz trace primarily the cold dust 
component, it in particular will be tracing a more strongly 
compressed average density. For our dust model, we use the 
same parametrized spiral arm modulation but with a con- 
trast of 7. This value, however, is degenerate with other 
parameters such as the dust arm width. 

Our adopted contrasts are therefore either uncon- 
strained a priori or roughly compatible with the loose the- 
oretical constraint from the shock scenario. We will show 
that we can reproduce the degree of polarisation in the data 
in the outer Galaxy with our current choices, but we note 
that this will not provide a unique solution. Additional data 
on the distribution of both the dust and the fields in the 
Galaxy will be needed to create an unambiguous model of 
the effects of compression on the ISM in the spiral arms. 

4.3 Cosmic ray arm contrast 

The situation for CREs is somewhat different from that 
of the magnetic fields and dust. Similarly to the thermal 
gas, CREs can be considered to be tied to field linefl, so 

^ Although magnetic field lines are a conceptual tool rather than 
a physical reality, it is always possible in ideal MHD to define field 
lines in such a way that their perpendicular velocity coincides with 
the flui d perp e ndicu lar velocity, itself equal to the electric drift 
velocity ISternI l ll966l) . Typical CREs, which also have a perpen- 



that they should undergo the same perpendicular compres- 
sion as the field. However, in the realistic situation where 
the shock front is not infinite and/or the field has a small 
component across the shock, CREs would diffuse along field 
lines at roughly the Alfven speed, thereby reducing their 
post-shock density. A n additional effect, clearly discussed in 
iFletcher et al.l (|201lh . is that the perpendicular momentum 
squared of each CRE should increase at the shock in direct 
proportion to the field strength - by virt ue of the conser - 
vation of the first adiabatic invariant (e.g.. lNorthrod[l963l ). 
This would lead to an overall shift of the CRE energy spec- 
trum towards higher energies, and hence to an additional 
increase (typically by the same factor of 4 as for the gas and 
the magnetic field) in the number density of CREs with a 
given energy, e.g. the characteristic energy for a given radio 
synchrotron observation. 

From these two different effects on the CREs, each of 
which may lead to a factor of 4 compression, the number 
density of CREs at a given energy is expected to increase 
at the spiral-arm shock by a factor between 4 and 16. This 
abrupt increase in density would then be smoothed out on 
the downstream side, on roughly a diffusion length, L ~ 
%/ Dt, where D is the diffusion coefficient (usually assumed 
to be isotropic for simplicity) and t is the CRE lifetime. 
With D ~ 3 X 10^* cm^ s"^ and t ~ (lO'^ - 10*) yrs for 
electrons of a few GeV in a field of a few \iG, the diffusion 
length works out to be L ~ (1 — 3) kpc. 

Another aspect that might affect our analysis is the 
question of anisotropic diffus i on of CRs from their acceler- 
ation sites. lEffenberger et al.l (|2012l ). for example, simulate 
the propagation of CR protons including anisotropic diffu- 
sion in the presence of a spiral magnetic field. They find 
that if the CR sources are also distributed along spiral arms, 
the resulting CR density retains a strong azimuthal modu- 
lation, with an amplitude ratio as large as a factor of ~ 6 
for strongly anisotropic diffusion (as opposed to only ~ 2 for 
isotropic diffusion) in the case of 1 GeV protons. Note that 
this effect is completely separate from, and not necessarily 
spatially coincident with, the CRE modulation due to the 
spiral-arm shock itself. 

The question of whether the CRE distribution varies 
from arm to inter-arm regions is therefore a relevant uncer- 
tain ty in our analysis. 

lAckermann et al] (|201ll ) present a study of the third 
Galactic quadrant where comparison of the 7-ray data with 
the velocity cubes in HI allows a separation of the spiral arms 
from the inter-arm regions. Among their results is a measure 
of the 7-ray emissivity in each region. They find that the 7- 
ray emission properties, when normalised by the varying gas 
density, vary little as a function of Galacto-centric radius. 
The small (« 5 per cent) average difference between the 
arm and inter-arm region implies a slightly lower emissivity 
per unit gas density in the Perseus arm compared to the 
inter-arm region that lies just inward of it. Though this is 
an analysis based on CR nucleon interactions rather than 
CR electrons, there is as yet no evidence in favour of any 
CR enhancement in the Perseus arm. 

Since no enhancement as large as predicted has so far 



dicular velocity equal to the electric drift velocity, can then be 
regarded as moving with the field lines. 
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been ob served either toward our Galaxy's Perseus arm or 
in M51 l|Fletcher et al.ll201lh . we consider our model of a 
smooth, azimuthally symmetric CRE distribution to be suf- 
ficiently realistic for this first-order analysis. We discuss the 
small effect a CRE enhancement in the arm might have on 
our results in § O 



4.4 Dust intrinsic polarisation fraction 

The details of the thermal dust emission process, and in par- 
ticular the alignment of the grains wi th the magnetic field, 
are not well understood jPraine fc Frai ssc 200!|). There are 
therefore few theoretical constraints on the intrinsic polari- 
sation fraction of dust, in contrast to the case of synchrotron, 
which introduces an additional uncertainty into the analysis. 
Observations of nearby, high-latitude clouds, which are likely 
observed at nearly their intrinsic polarisation, have dust po- 
larisat ion fractions of up to ~20 per cent with most much 
lower (jBenoit et al.ll2004l : [Ponthieu et al]|2005l : lFauvet et all 

[ml. 

As we will discuss further below, the interesting differ- 
ences between our model and the observations include the 
overall level of dust polarisation. We will show that we can 
address this problem in the outer Galaxy by changing the 
spatial distribution of the magnetic field components rela- 
tive to each other and to the dust emitting regions. In other 
words, a higher than expected polarisation fraction can be 
reproduced by increasing either the intrinsic polarisation or 
the degree of ordering in the fields in the dust emitting re- 
gions; the two questions are degenerate. Adopting a higher 
degree of intrinsic polarisation is therefore the conservative 
choice that avoids over-stating the data's ability to constrain 
the field ordering. 

As described in § 13.21 our simulation has a finite resolu- 
tion that in this large-scale study is of order 20 pc. The 
number we use as the "intrinsic" polarisation fraction is 
therefore the assumed average degree of polarization from 
a region of this size. Since the average molecular cloud size 
is of the same scale, this fraction then represents an average 
over a cloud. We might improve the modeling of the sub-grid 
effects, i.e. include the turbulence driven within the clouds 
with a model that could be significantly different in both 
scale and amplitude from what we have modeled for syn- 
chrotron. This improvement, however, could only decrease 
the predicted degree of polarisation from that obtained as- 
suming an intrinsic polarisation of 30 per cent. 

We therefore adopt 30 per cent for the intrinsic polar- 
ization fraction as a conservative upper limit while noting 
that this value is degenerate with the question of how or- 
dered the fields are in the dust emitting regions. 



4.5 External galaxies 

External galaxies have been mapped in both total an d po- 
larised intensity by low-frequency radio observations. iBeckl 
l|2009D reviews the different morphologies seen and in partic- 
ular, the different relationships observed between the spiral 
structure seen in polarised synchrotron emis sion and tha t 
seen in gas tracers. A study of M51 bv .Fletcher et al.l (|201ll ) 
provides one of the first estimates of arm versus inter-arm 
contrasts seen in polarised synchrotron emission. They find 



a much lower level of contrast than expected by the scenario 
of shock compression by the spiral wave, but they note that 
it is unclear whether their observations have sufficient reso- 
lution. Their limit on the size of the compressed region based 
on the observed low level of contrast and the telescope beam 
is n ot so small as to rule o ut the scenario. 

iPatrikeev et all l|2006l ) use a wavelet analysis to find 
the arm ridges in each of a set of maps from low-frequency 
radio polarisation to CO and infra-red observations. They 
find a systematic shift in the molecular gas versus the mid- 
infrared in each spiral arm, consistent with the delay of star 
formation triggered in the shock. The polarised emission is 
also the component the furthest upstream, as expected if the 
magnetic fields are compressed and aligned with the shock. 

Though these observations do not provide strong con- 
straints on the relative positions or widths of the arms or 
on the arm versus inter-arm contrasts, they do provide evi- 
dence for such shifts and for variations between components' 
arm geometries. As such, observations of external galaxies 
motivate our attempt to explore these differences in our own 
Galaxy. 



5 RESULTS 

5.1 Profiles along the Galactic plane 

Following IJaffe et all l|2010l ). we examine profiles of the ob- 
servables along the Galactic plane as shown in Fig. [T] The 
model in b lue is simp l y an updated form of our previous 
model from lJaffe et all |201ll ). The RMs constrain the coher- 
ent field intensity in each arm, while the total and polarised 
synchrotron intensities constrain the amounts of isotropic 
and ordered random components. 

We have added the dust emission as described in § 13.11 
and see that the model (blue curve in Fig. [l| does not fit 
the dust polarised emission data at all, with two main prob- 
lems. Firstly, the shape of the profile is not correct, with 
proportionately too much polarised emission in the inner 
Galaxy compared to the outer Galaxy. Secondly, in the outer 
Galaxy, the polarisation fraction is too low. 

It is apparent from t he RM panel of Fig. [1] that the 
data from Ivan Eck et all l|201ll ) do not match our model. 
The geometry of our model was originally constrained in 
Ijaffe et al.l l|2010l ) before these data were available. It is clear 
that though the amplitude of the closest peak in our RM 
model roughly matches the van Eck et al. data, its longitude 
does not. We note that simply changing the orientations of 
the spiral arms, their pitch angles (currently —11.5°), or the 
radius of the inner Galaxy ring does not solve the problem. 
Something more complicated than simple logarithmic spirals 
and an annular ring will be needed, but the average RM 
and its large scale variations predicted by our model are 
sufficiently accurate to fix the strength of the coherent field 
component. 

It must also be noted that there is significant uncer- 
tainty regarding the distance to the Perseus arm in the outer 
Galaxy. The spiral arm model used in the NE2001 model 
places it roughly 2-3 kpc away toward the anti-c entre, while 
the analysi s of velocity data presen ted in both lAbd o et al.l 
(2010) and lAckermann et al.l |20l3) imply nearly twice the 
distance. There is a large degree of uncertainty in both esti- 
mates, but our analysis is not significantly affected by this. 
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Our results are dependent only on the relative positions of 
the diflFerent components, not their absolute distances, and 
can therefore be re-scaled to match a corrected distance de- 
termination. 



5.2 Outer Galaxy 

For these simulations, we have assumed an intrinsic polar- 
isation fraction of the dust of 30 per cent. As discussed in 
Section |4]4l 30 per cent is at the high end of what is consid- 
ered plausible and significantly higher than what is usually 
assumed. By using 30 per cent, we are therefore being con- 
servative, since a more commonly assumed 15 per cent makes 
the discrepancy between the model and the data in the outer 
Galaxy more extreme. An even higher intrinsic polarisation 
fraction might solve the problem in the outer Galaxy but 
would be in even more serious confiict with observations. 

Nor can we reproduce the observed polarisation fraction 
by increasing the average amount of magnetic field ordering, 
because this would break the fit to the synchrotron emission. 
But since the synchrotron emission depends on the spatial 
distribution of CREs, and since these CREs quickly diffuse 
away from their acceleration sites (see discussion in Sec- 
tion I4.3|l , the synchrotron polarisation degree reflects only 
the average level of coherent and ordered random fields ver- 
sus isotropic random fields. It is not very sensitive to where 
those components are spatially located. Therefore, we can 
shift the components relative to one another without much 
changing the synchrotron emission profiles. 

We show the effect of this with a new model in yellow 
in Fig. [2l which is identical to the previous model except 
for the location of the isotropic random component. In the 
blue curve, as before, all field components as well as the 
dust have the peak of the arm in the same place. But the 
dashed yellow curve is the prediction from a model where the 
coherent and ordered random fields peak in the same place 
as the dust but the isotropic random component has been 
shifted away by 45° in Galacto-centric azimuth. This places 
it exactly between two dust arms, the furthest separation 
possible, maximising the resulting polarisation. The model 
is shown in Fig. |3] The result is two very similar profiles for 
the synchrotron polarisation fraction but significantly more 
polarised dust emission, since it is then emitted in regions 
of strongly ordered fields. This then matches far better the 
data in the outer Galaxy. 

The model parameters have been set to those that best 
match the outer Galaxy, particularly the third quadrant 
(longitudes 180° to 270°), since it is unclear how much the 
so-called "Fan" region (dominating most of the second quad- 
rant, 90° to 180°), which remains anomalous, may be due 
to a local feature. 



5.3 Inner Galaxy 

The shift of the isotropic random field component away 
from the dust emitting regions makes the mismatch between 
model and observations in the inner Galaxy worse. The in- 
ner Galaxy is also complicated by the presence of what is 
often referred to as the "molecular ring" and whose effect on 
the magnetic field components is not known . Furthermore , 
our Galaxy is thought to have a bar (see, e.g.. lGerharclll2003 ) 



with again unknown effect on the field components and on 
the dust. 

The parameters that match the dust emission in the 
outer Galaxy fail in the inner Galaxy, and this tension can- 
not be resolved with the current parametrization of the 
fields. Clearly, a more complicated model that smoothly 
varies quantities such as the arm positions, widths, con- 
trasts, and pitch angles from the inner to outer Galaxy will 
be necessary to fit the data. 

5.4 Comparison to low-frequency radio 
polarisation 

Note that the model in Fig. |3]has the Sun position (8.5 kpc 
to the left of centre) just inside the isotropic random field 
arm. This means that looking toward the outer Galaxy, 
we ought to see Faraday depolarisation effects in the low- 
frequency radio data in that direction. Though the strength 
of the depolarisation effect is unknown and would require 
some additional modeling to estimate, it is clear that the 
effect is the opposite of what is seen, namely that it is the 
inner Galaxy where there is significant Faraday depolarisa- 
tion and the outer Galaxy o nly where strong polar isation 
is apparent. See, e.g., fig. 1 of iBurigana et al.l (|2006l ) show- 
ing the polarised emission at 1.4 GHz. The implication is 
that the positions of the arms in our model may need to be 
adjusted relative to the position of the Sun. 

We will in future add to our studies the Faraday depo- 
larisation effects seen in the 1.4 GHz polarisation map, since 
this is important additional data on the spatial distribution 
of the turbulent medium. 



6 DISCUSSION 

We show in Section [5] how the profiles of polarised syn- 
chrotron and dust emission in the Galactic plane, combined 
with previous work to constrain average strengths of the field 
components in the Galaxy, require a spatial separation of the 
different field components. It is clear from Fig.[2]that a con- 
stant mixing of the coherent plus ordered random fields with 
the isotropic random component is inconsistent with the po- 
larised dust data in the Galactic plane. We can also see that 
the problem is opposite in the inner and outer Galaxy. Our 
current model cannot simultaneously fit both regions, but 
we can already draw some conclusions from Fig. [2] 

6.1 The shock scenario 

Firstly, we can reproduce the emission in the outer Galaxy 
with a model where the isotropic random component peaks 
in arms that are clearly separated from both the coherent 
and ordered random components and the dust emitting re- 
gions. The synchrotron emission constrains only the average 
amount of coherent plus ordered random and isotropic ran- 
dom magnetic fields in the Galaxy, but does not strongly 
depend on where they are located. External galaxies show 
indications of magnetic arms distinct from the gas arms and 
with morphologies th at are not nec essarily the same (see, 
e.g., iFrick et al.l |2000| or iBeckl [20091 ') . They can be shifted 
relative to each other, have different pitch angles, or cross 
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Figure 2. Panels as in Fig.[T] The blue line is the model from Fig.[T]and has coincident arms for all components while the dashed yellow 
model is that where the isotropic random field component has been shifted back 45° in Galactic azimuth, as shown in Fig. \3\ See § \5\ 



over the gas arm ridges. There is clear evidenc e for such ef- 
fects in the case of M51 jPatrikeev et al.ll2006l ). Our results 
in the outer Galaxy support the idea of a shift between com- 
ponents. 

The second conclusion we can draw from Fig. [2] is that 
the situation is reversed in the inner Galaxy, where the data 
show no peak in the polarised emission as predicted, im- 
plying that the fields are less ordered in the dust emitting 
regions than in our models. Our current model cannot be 
tuned to simultaneously fit both regions. 



We can reproduce the polarisation in the outer Galaxy 
with a model that has the peak of dust emission spatially 
coincident with the coherent and ordered random fields but 
offset from the isotropic random field. Such a lag can be 
naturally explained by the idea of a large-scale shock wave 
associated with the spiral arm: the shock compresses the 
fields so they lie predominantly parallel to the shock, com- 
presses the gas and dust likewise, and triggers collapse of 
dust clouds and subsequently star formation. The turbu- 
lence is then generated in supernovae downstream of the 
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Figure 3. Top panel: contour plot of the Galaxy model seen from above. The magenta traces the coherent and ordered random magnetic 
fields. In blue is the dust distribution. In gold is the isotropic random field component. The Sun position is at -8.5 kpc, to the left of the 
galactic centre. The lower panel shows the profile on the axis connecting the Sun and the Galactic centre. The observed profiles along 
the plane are shown in Fig. [2] 



compression after the stars have had time to evolve, tlius 
explaining the shift in the isotropic random field component 
relative to the other arms. 

This scenar i o is also c onsistent with the studies by 
iHaverkorn et all (|2006l . I 20081 ) . which examine the structure 
function of Faraday RM fluctuations in directions either tan- 
gent to spiral arms or tangent to inter-arm regions. They find 
that in the inter-arm regions, the structure function shows 
a rising spectrum up to a cut-off scale, consistent with tur- 
bulence from a Kolmogorov cascade of energies from large 
scale turbulence injected by supernova remnants. The arm 
tangents by contrast do not show such a spectrum but rather 
a flat structure function, perhaps implying a much smaller 
outer scale of turbulence, small enough that the turn-over on 



the structure function is not visible. This is interpreted by 
Haverkorn et al. as an alternative source of turbulence, e.g. 
HIT regions and stellar winds instead of supernova remnants. 
But it is also what one might expect to see if the larger scale 
turbulence has been compressed perpendicular to the shock, 
i.e. perpendicular to the LOS when looking tangent to the 
arms. Alternatively, it could simply be that the the star for- 
mation is triggered in the shock and only later evolves into 
the supernova remnants that inject the turbulence, which 
would mean the turbulence is always downstream of the 
compressed region. 

Since the downstream region is leading the shock inside 
co-rotation (where the interstellar gas rotates faster than 
the spiral pattern) and lagging behind the shock outside co- 
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rotation, our simple model of spiral arms being azimuthally 
shifted as wholes is not accurat e. The value of t he co-rotation 
radius is very uncertain (e.g., lGerhardll201ll ). but it could 
be close to the Galacto-centric radius of the Sun, in which 
case inner/outer Galaxy would roughly coincide with in- 
side/outside co-rotation. In other words, the depolarising 
turbulence would be shifted ahead of the dust arms in the 
inner Galaxy, and behind the dust arms in the outer Galaxy. 
But this difference in the shift direction would not change 
the polarisation degree. 

It should be noted that the dust is emitted mostly in 
relatively compact clouds that are not resolved by our simu- 
lations. One can easily imagine different degrees of ordering 
in the magnetic fields in the clouds compared to the dif- 
fuse ISM traced by synchrotron emission, e.g. by turbulence 
injected at smaller scales. As discussed in § 14.41 however, 
our modeling effectively assumes 30 per cent polarization 
at the scale of a typical cloud, which is significantly higher 
than anything observed. Furthermore, there is also a truly 
diffuse dust component that would trace the large-scale tur- 
bulence seen in synchrotron emission, though the ratio of 
such emission to that from relatively compact clouds is un- 
known. We therefore find it difficult to envision a plausible 
alternative scenario to increase the dust polarization in the 
outer Galaxy as long as the dust clouds remain embedded in 
the same turbulent ISM that is traced by the synchrotron. 

6.2 Refinements for the inner Galaxy 

There are several possibilities that might explain the differ- 
ence in the inner Galaxy. Broadly speaking, they are either 
that the dust-emitting regions have more disordered fields in 
the inner Gala:xy, or that the intrinsic polarisation fraction 
of the dust emission drops. 

We could perhaps model the former by allowing the 
field components to have different pitch angles so that the 
isotropic random component is coincident with the dust in 
the inner Galaxy but gradually shifts away in the outer 
Galaxy. This would result in a degree of ordering that would 
vary with Galacto-centric radius and might be made to re- 
produce the data. But this would not be consistent with the 
shock compression scenario that naturally accounts for the 
polarisation in the outer Galaxy, and we currently lack a 
plausible alternative scenario for such a morphology. 

There is similarly a question of when the collapsing dust 
clouds decouple from the magnetic fields. There is more star 
formation in the inner Galaxy, and the colder dusty regions 
could be less strongly coupled to the large-scale ordered 
fields in the inner Galaxy. Even if the clouds remain coupled 
to the fields while collapsing, the collapse itself could disrupt 
the ordering in the fields. Since the synchrotron emission 
comes from a large volume while the dust emission may be 
dominated by relatively isolated regions of collapsing clouds, 
the dust emitting regions may have more random fields at 
small scales embedded in the larger-scale fields traced by 
synchrotron. 

Alternatively, the truly intrinsic polarisation fraction 
could drop in the inner Galaxy due to the different environ- 
ment (e.g. temperature) of the dusty regions. As discussed 
in § 14.41 the intrinsic polarisation fraction is not well under- 
stood. 

We could fold these possibilities into an effective sub- 



grid polarisation fraction that varies with Galacto-centric 
radius. Our large-scale analysis cannot distinguish these pos- 
sibilities without additional external constraints on small- 
scale dust polarization properties, but higher-resolution 
studies of the fields in individual clouds may be very in- 
formative. 

6.3 Caveats 

Note that the uncertainty in the true intrinsic polarisation 
fraction of dust does not affect our analysis in the outer 
Galaxy. As discussed in § 14.41 it is very unlikely to be high 
enough to account for the observed polarisation without spa- 
tially separating the dust from the isotropic random fields. 

Another uncertainty in this analysis is the impact of 
the possible modulation of the CRE density due to the two 
effects discussed in § 14.31 namely the compression and ac- 
celeration of the CREs at the spiral-arm shocks and the 
distribution of CRE injection sites. Both sources of CRE 
density modulation arise in spiral arms, but they are not 
necessarily spatially coincident, as CREs are presumably in- 
jected by supernova remnants in the turbulent regions that 
develop downstream from the shocks, away from the ampli- 
fied fields. One can speculate that the two effects may result 
in a fairly uniform average CRE density across spiral arms, 
thus explaining the observed lack of any such modulation in 
the CRE density. 

In order to reproduce the observed level of synchrotron 
emission, and not to result in arm tangents that are more 
peaked in emission than observed, any CRE enhancement 
would have to be offset by a corresponding lowering of the 
coherent and ordered random magnetic fields. In any case, 
the degree of polarisation of the dust emission does not de- 
pend on the strength of the fields, only on the degree of 
ordering in the dusty regions. So the question of the CRE 
density is an interesting one but does not have a large effect 
on our main conclusions. 

Neither the uncertainty in the intrinsic polarisation 
fraction of dust emission nor in the CRE density modula- 
tion, therefore, will change our main conclusion: that in the 
outer Galaxy, the dust emitting regions have more ordered 
fields than the average as traced by synchrotron emission. 



7 CONCLUSIONS 

We have continued our project to model the large scale 
Galactic magnetic fields by adding the polarised thermal 
dust emission from the WMAP 94 GHz band. We find that 
our model previously fit to synchrotron and Faraday RMs 
does not reproduce the observed level of dust polarisation. 
Even with an assumed intrinsic polarisation fraction of 30 
per cent, which is higher than observed, the level of predicted 
polarisation is too low in all but the innermost regions of the 
Galactic plane. This implies that the dust emitting regions 
have more ordered magnetic fields than average. We are able 
to reproduce the data in the outer Galaxy by modifying the 
model to spatially separate the isotropic random field com- 
ponent from the dust emitting regions and the coherent and 
ordered random field components. This separation is consis- 
tent with the scenario of fields and dust compressed in the 
spiral arm shocks, with turbulence generated by supernova 



© RAS. MNRAS OOO.ITlfTSl 



12 T. R. Jaffe et al. 



remnants downstream. The same model results in an over- 
prediction of the dust polarisation in the inner regions of the 
Galaxy, which indicates that a more sophisticated model is 
necessary. We will in further work explore models that allow 
either the spiral arm characteristics of each component or 
the dust intrinsic polarisation degree to vary with Galacto- 
centric radius. Our models would be improved by numer- 
ical simulations of galaxies' magnetic fields that included 
arm versus inter-arm variations and provided constraints to 
break some of our parameter degeneracies. This project will 
be significantly aided by the Planck satellite's more sensi- 
tive observations of polarised dust emission at both large 
and small scales. 
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APPENDIX A: MODEL PARAMETRIZATION 

The following tables describe the parametrization of the 
magnetic field and dust distribution models. The difference 
to previous works is in simplifying the radial profiles (since 
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we do not fit t he in nermost region of the Galaxy where 
the iBroadbentI l|l989l ) features have effect) and in making 
the field continuous in the transitions between arms. Each 
component has an associated scale height, which is uncon- 
strained by our analysis confined to the Galactic plane. But 
they are given for completeness and for the small effect they 
may have within the resolution of 6° . 
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Total field Btot = Bcoh + Biso + Bord 
Param. Default Equation Description 



Axisymmetric spiral coherent magnetic field: BassC^": </'i ^) = S(r)_Bcoh (2)6 



Bo 


1 n-G 






Global amplitude normalisation 


^scalc 


20 kpe 


i?(r) = Bocxp(-rVi?4^lJ 




Outer radial profile parameter. See fig. 3 of 
1 Jaffe et al.l ll2010f). 




6 kpe 


Bcoh{z) = sech^(z/hg) 




Global scale height of the coherent field. 


^max 


20 kpe 






Maximum radius, beyond which |B| = 


^mol 


5 kpe 


g_r (eos(0-|-7r/2)e:,,sin(</.-|-7r/2)ey,O) if 
\ (sin(ep )er , eos (9p)e^ , 0) if 




Field direction circular within "molecular 


r > Rniol 


ring" and otherwise spiral (expressed in cylin- 
drical coordinates, with the x-axis between the 
Sun and Galactic centre). 


Op 


-11.5° 






Op is the pitch angle of field direction and spiral 
arms 






Spiral compression arm coherent magnetic field: Bcoh = 


= Bass + l^f^"™" Barm.i 


di 


(3,0.5, -4, 
1.2,-0.8) 


Barm.i = B{r)ai pc{di)'B 

pc{d) = c{r)Bcomp{z) cxp{-{d/do{r))'^) 




Amplitude of each of four spiral arms and 
molecular ring.. 

Additional compression factor relative to back- 
ground, d is the distance to the nearest arm in 
kpe, computed using ri((j>) 


4'o,i 


10° +90°i 


r,{<l>) = Rs exp [((/) - 0o,i)//3] and /3 = l/tan{6ip) 




r{cp) gives arm radius at given azimuth, and 
00,1 the azimuthal orientation of the rotation 
of spiral around axis through Galactic poles 


-^^comp 


7.1 kpe 






Scale radius of compressed spiral arms. 


Co 


2.5 


. N / Co if r ^ rcc 
\ Co(r/rcc)-^ ifr>rcc 




Arm compression amplitude, tailing off after 
a region of constant compression. The con- 
trast, e = n^'*) /n^^\ is the ratio of downstream 
(post-shock) to upstream (pre-shoek) density, 
in which ease Co ~ e — 1 if we consider arm 
versus inter-arm instead of down- versus up- 
stream.) 


^cc 


12 kpe 


See c(r). 




Region of constant compression. 


do 


0.3 kpe 


do(r) = do/(c(r)B(r)) 




Defines the base width of arm enhancement, 
which varies with radius. 


he 


2 kpe 


Bcomp(z) = sech^(z//ic) 




Scale height of the spiral compression. 






Isotropic random and ordered random magnetic fields: 






Bi,o + Bo,d = [BgrfBgrf(2)Bgrf{'") + 


^—-i arm. 


1 y^iV -Dord 
i ^—Ji arm.i 


h^rns 


2 kpe 


BgRf{z) = SCcil^iz/hirns) 




Global scale height of the turbulent field field. 




20 kpe 


BGRFir) = exp(-r^/Rl.j^p) 




Global scale radius of the turbulent field. 


a 


-0.37 


Peik) = (|BGRF(fc)P) <x fc° and 

BGRF(x)=J^-MBGRF(k)] 




Power law spectral index of initial GRF; de- 
fault of -0.37 from iHan et al.l \2004) in ID 
(compare to 3D Kolmogorov turbulence pre- 
diction of —5/3). 


Deo 


0.1 kpe 


BGRF(fc) = for fc < l/Dco 




Cutoff maximum of GRF fiuetuations (mini- 






Brms = (Bgrf(x)) ''' 




mum determined by resolution) 


Brms 


3.5 




Total RMS amplitude of GRF fluctuations 


ford 


0.15 


BLon.i = P-Bgrf and B-f^ . = /orfPcBp^.j 




Relates ordered to isotropic random compo- 
nent. Bproj is the component of the GRF par- 



allel to the coherent field direction, i.e. the spi- 
ral. 



Table Al. Table of magnetic field modelling parameters as described in S 13.11 The total field is the sum of the coherent, isotropic 
random, and ordered random components. The coherent and isotropic random each have a (statistically) axisymmetric and a spiral 
component, while the ordered random only has the spiral component. Parameter values for the spiral are the same for each component 
in this work with the exception of Cq, which in Fig. [2] is different for the isotropic random component only. 
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Dust emissivity model parameters: 


Parameter 


Default 


Description 


her A 


0.4 kpc 


Scale height of dust uncompressed component. 


hc,d 


0.1 kpc 


Scale height of the dust arms. 


^scale,d 


8 kpc 


Scale radius of the uncompressed component. 


^comp,d 


13 kpc 


Scale radius of the compressed spiral arms. 


Co,d 


6 


Compression factor for dust. 


1^0, d 


0.1 kpc 


Base width of dust arms. The width varies with 




r as described in Table \KT\ 



Table A2. Table of parameters for the dust density distribution. 
The parametrization is the same as for the magnetic field models 
defined in Table lAll an axisymmetric component plus a spiral 
compressed component, and any parameters not listed here are 
the same. 
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